Decoherence of two qubits in a non-Markovian squeezed reservoir 
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The decoherence of two initially entangled qubits in a squeezed vacuum cavity has been investi- 
gated exactly. The results show that, first, in principle, the disentanglement time decreases with the 
increasing of squeeze parameter r, due to the augmenting of average photon number of every mode 
in squeezed vacuum cavity. Second, there are entanglement revivals after complete disentanglement 
for the case of even parity initial Bell state, while there are entanglement decrease and revival before 
complete disentanglement for the case of odd parity initial Bell state. The results are quite different 
from that of the case for qubits in a vacuum cavity. 
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In recent years, the entanglement dynamics of qubits, 
coupled with environment, has attracted much atten- 
tion. Many works have been devoted to the phenomenon 
termed as "entanglement sudden death" (ESD)[l|,|3,[l,|3, 
S i 0, S 1 mM, El El, Q [II- it is shown that 
spontaneous disentanglement, for two-level atom model, 
may take only a finite-time to be completed, while lo- 
cal decoherence (the normal single-atom transverse and 
longitudinal decay) takes an infinite time[3]. And, for a 
non-Markovian reservoir of initially vacuum cavity, there 
is revival phenomenon, which allows the two-qubit entan- 
glement to reappear after a dark period of time, during 
which the concurrence is zero|12l|. In general, the char- 
acteristic of the environment plays an important role in 
the evolution of multi-particle entanglement. Up to now, 
many typical environments have been investigated, such 
as, vacuum, squeezed vacuum, multimodc vacuum cavity, 
single mode cavity and so on. 



Recently, cavity systems with very strong couplings 
have been discussed 16]. Generally, in atom-field cavity 
systems, this ratio is typically of the order 10~^ ~ 10~^. 
However, the ratio may become order of magnitudes 
larger in solid state systems, and the full Hamiltonian, 
including the virtual processes (counter- rotating terms), 
must be considered [17|. In this paper, we will focus on 
the decoherence of two qubits strongly coupled with a 
non-markovian squeezed vacuum reservoir by the method 
in Ref . p^ . In section. 2, the reduced non-perturbative 
non-Markovian quantum master equation of atom is de- 
rived and its exact solution is obtained. In section. 3, 
the decoherence of two initially entangled atoms, coupled 
with two cavities separately, has been discussed. The 
conclusion is given in section. 4. 



A. Hamiltonian and non-perturbative master 
equation 

Now we restrict our attention to two noninteracting 
two-level atoms A and B coupled individually to the cav- 
ity field. To this aim, we first consider the Hamiltonian 
of the subsystem of a single qubit coupled to its reservoir 
as 



where 
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where ujq is the atomic transition frequency between the 
ground state |0) and excited state |1). Cz = |1)(1|~|0)(0|, 
(T+ = |1)(0| and cr_ = |0)(1| are pseudo-spin oper- 
ators of atom. The index k labels the field modes 
of the reservoir with frequency Wfc, a], and Ofc are the 
modes' creation and annihilation operators, and gk is 
the frequency-dependent coupling constant between the 
transition |1) — |0) and the field mode k. 

The reduced non-perturbative non-Markovian quan- 
tum master equation of atom could be obtained by path 
integral El 



d /■* 

-j:Pa = -'iS,aPa ~ / ds{2.a 
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where Zq, Za and Zar are Liouvillian operators defined 
as 
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£oP = 


[Ha+Hr,p] 


ZaP = 


[Ha,p] 


^arP = 


[Har,p] 



and (...)r stands for partial trace of the reservoir. 

Then, we assume the reservoir is initially in squeezed 
vacuum state (l9ll 
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Sk = exp{re '■'^ak^+kak.-k - re'^al^+k'^l.-k) C^) 

r, 9 are squeeze parameters, |Ofc) is the vacuum state of 
mode k. The central frequency of squeezing device is 
Wo = ckc, which corresponds to multimode squeezed vac- 
uum state with the central frequency equal to the cav- 
ity resonance frequency and atom transition frequency. 
And the spectral density of the reservoir is in Lorentzian 
formlH, m 



v+ = 2[M"'f + M"a^ + M^a'+Nf+{N+l)a^] 
V- = 2[M"'f + M'^a^ + M^a^ + {N + l)f + Na^] 
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-[1 - exp{--ft)] 



N = sinh'^r, M = -e"^ sinhrcoshr. a^, a^ , a* and M^, 
M^ , M* are real part, image part and conjugate of a and 
M, respectively. Here, N is the average photon number 
of every mode in squeezed vacuum cavity. And M, M* 
represent the phase-dependent correlation between dif- 
ferent modes as {hkhk') and {b\b\,), resDectivelv|19j. 



J(u;) 



gl{5{uj - ujk) + 5{lo + LOk)} 
1 A72 



k 



(8) 



27r {uj - t^o)^ + 7^ 

where 7 represents the width of the spectral distribution 
of the reservoir modes and is related to the correlation 
time of the noise induced by the reservoir, r^ = I/7. 
The parameter A is related to the subsystem-reservoir 
coupling strength. There are two correlation functions in 
this modelim 
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ai comes from the rotating-wave interaction and 0:2 from 
the counter-rotating wave interaction. 

So, the non-perturbative master equation of the sub- 
system could be derived from Eq.® 



d 



dt 



Pa 



-Tpa + [SqJo + e+J+ + e-J-]pa 

+ [voK^ + V+K+ + v_KJ\pa (11) 

where Jo, J+, J-, A'o, K+ and K- are superoperators 
defined as 

{(J+a_Pa + Pacr+a. 
<7+Pa<y~ 
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and 



r = 2M'^/ + 2(A/^a^ + AfV) + (2iV + l) (/ + «-") 

e+ = 2Mf + 2M*a+{2N +l){f + a) 
e- = 2Ma* + 2M*f + {2N+l){f + a*) 



B. Exact solution of master equation 

The time evolution of density operator in Eq. lllip could 
be obtained with algebraic approach in Ref.[22] because 
the superoperators herein satisfy SU{2) Lie algebraic 
communication relations, i.e. 

[J-,J+\pa = ~2JoPa 

[Jo,J±]Pa = ±J±Pa 

[K^,K+]pa = ~2KoPa 

[Ko,K±]pa = ±K±Pa 

[K,,J,] = (12) 

where i,j = 0, ±. By directly integrating Ea. (fTT|) . the 
formal solution is obtained as 12311 



p^(t) — g-^krpgfo dt{eoJo+e+.I++e-J-) 

where T is time ordering operator and 

Tk = (2A/^ + 2A^+l)(F + 5^) + 2A/^5^) 
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where a^ , d^ and d* are real part, image part and con- 
jugate of a, respectively. 

The exponential functions of superoperators in Eq. (|13| 
could be disentangled in form: 23] 



J'g/o dt{eoJo+£+J++£-J-) 
fglo dt{voKo+v+K++u-K-) 
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where j+, jo, j- and fc_(-, fco, fc_ satisfy the following 
differential equation 



X+ ^ p+- p-X_^ 
X(3 = Pq- 2/i_X+ 
X- — p-.exp{XQ) 
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/i = e for X = j and ji = v for X = k. 

Using the results above, the exact solution of the mas- 
ter equation Eq. (fTT|) is obtained 
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y - e-J«/2j+ (24) 
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FIG. 1: Concurrence C<s> as a function of ^t and 6 with A 
IO7, ujo = IO7, /3 = 72/2 and r = 0.2. 



In order to investigate the entanglement dynamics of 
the bipartite system, we use Wootters concurrence [2^. 
For simplicity, we assume that the two subsystems have 
the same parameters. The concurrence of the whole sys- 
tem could be obtained 



ci — 2e 
C2 = 2e 



:{0,Ci,C2},(^ = $,*) 
{\/P23P32 - \/pilPii) 



(25) 



-2r, 



-2r, 



{y/PliPil - VP22P33) 



corresponding to the initial states of |$) — /9|01) + ?7|10) 
and l^f*) = /3|00) + ?7|11), respectively. Where (3 is real 
and < ^ < 1, r; = \r]\e''^ and (3^ + \r]\^ = 1. For 
maximum entanglement Bell state, /3 is equal to V2/2. 
The reduced joint density matrix of the two atoms, in the 
standard product basis <B = {|1) = |11), |2) = |10), |3) = 
|01), |4) = |00)}, could be obtained by the method in 
ref.nj 



Pll pi4 

„AB _ „-2rt I P22 P23 
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here the diagonal elements are 

Pii = l'^pii{0) + lmp22{0) + mlp33{0) + m'^P44{0) 

P22 == IppiiiO) + lmp22{0) + mpp33,{0) + mnp44{0) 

P33 = lppii{^) +pmp22{0) + nlp33{0) + nmp44{0) 

piA = p^pii{0) +pnp22{0) + npp33{0) + n^p44(0) 

and the nondiagonal elements are 

/9l4 = x^PiiiO) +xyp23{0) + yxp32{0) +y'^P4i{0) 

P23 = xrpi4{0) + xqp23{0) + yrp32{0) + yqp4iiO) 

P32 = rxpi4{0) + ryp23{0) + qxp32{0) + qyP4i{0) 

P41 = r^Pi4{0) + rqp23{0) + qrp32{0) + q^P4i{0) 
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FIG. 2: Concurrence C* as a function of jt and /3^ with 
A = IO7, UJO = IO7, P = v/2/2 and r = 0.2. 



III. 



NUMERICAL RESULTS AND DISCUSSION 



In order to study the effects of non-Markovian squeezed 
reservoir on the decoherence, we assume that A is princi- 
pally equal to IO7 in Eq.®, which could be realized in 
a high-Q cavity l^j. 

First, we focus on the effects of squeeze parameters on 
the decoherence of two qubits with initial maximum Bell 
entanglement states |$) and |\I'), respectively. 

(A) For Wo = IO7 and /3 = v^/2, Fig.l and Fig.2 show 
that, for the case of initial maximum Bell state |$), the 
concurrence first decreases to a certain value and then re- 
vives before it vanishes, while that periodically vanishes 
and revives with a damping of its revival amplitude for 
the case of initial maximum Bell state | ^) . It also reveals 
that the amplitude and the duration time of entangle- 
ment revival are different for the case of tt/2 < < 3tt/2 
and for the case of < 6* < 7r/2 and 37r/2 <9<2tt. 

(B)For LOo = IO7 and /3 = ^2/2, Fig.3 and Fig.4 show 
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FIG. 3: Concurrence C* as a function of 'yt and r with A 
IO7, uJo = IO7, 13 = V2/2 and 6 = 7r/4. 



FIG. 5; Concurrence C* as a function of 7i and 13^ with 
A = IO7, u)o = 127, 61 = 7r/4 and r = 0.2. 





FIG. 4; Concurrence C* as a function of 7^ and r with A 
IO7, uio = IO7, /3 = V2/2 and 61 = 7r/4. 



FIG. 6: Concurrence C* as a function of 7i and (3 with 
A = IO7, uo = IO7, 61 = 7r/4 and r = 0.2. 



that, for r < 1, the concurrence first decreases to a cer- 
tain value and then revives before it vanishes for the case 
of initial maximum Bell state 1$), while that periodically 
vanishes and revives with damping amplitude for the case 
of initial maximum Bell state \^). And the concurrence 
decreases monotonically and vanishes permanently after 
a short time for r > 1. The results exhibit that the disen- 
tanglement time decreases with the increasing of squeeze 
parameter r for fixed 0, due to the increasing of average 
photon number of every mode in squeezed vacuum cavity. 

The above results also reveal that the decoherence be- 
havior of concurrence is sensitive to squeeze parameters 
and insensitive to initial state. 

Then, for fixed squeeze parameters 6 ~ 7r/4 and 
r = 0.2, the decoherence for different initial states was 
discussed. 

(A) The decoherence behavior of C3, is discussed as 
follows. 

(1) From Fig. 5, Fig. 6 and Fig. 7, we could find that 
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FIG. 7; Concurrence C* as a function of ^t and 0^ with 
A = IO7, u)o = 6.57, e = 7r/4 and r = 0.2. 
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FIG. 8: Concurrence C<s as a function of jt and l3 
A = IO7, Wo = 37, e = 7r/4 and r = 0.2. 



with FIG. 10; Concurrence C* as a function of yt and f3^ 
X = IO7, cjo = IO7, e = 7r/4 and r = 0.2. 
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FIG. 9: Concurrence C* as a function of jt and /3^ with 
A = IO7, cjo = 127, e = 7r/4 and r = 0.2. 



FIG. 11; Concurrence C* as a function of -yt and P^ 
A = IO7, Luo = 6.57, 61 = 7r/4 and r = 0.2. 



with 



the concurrence C$ will decreases and revives before it 
vanishes permanently. With the ratio of the coupling 
strength to the atomic frequency increasing, the times of 
entanglement revival decreases because of the increasing 
of the average value of correlation function 0:2, corre- 
sponding to the counter-rotating wave terms. 

(2) For Wo = 37 and A = 2O7, Fig.8 exhibits that 
the concurrence C$ decreases monotonically and vanishes 
permanently in a short time even in a non-Markovian 
squeezed reservoir, resulted from the strong interaction 
between atom and the non-Markovian squeezed reservoir. 

(B)The decoherence behavior of C3. is discussed as fol- 
lows. 

(1) From Fig. 9 , Fig. 10 and Fig. 11, we find that 
the concurrence periodically vanishes and revives with 
a damping of its revival amplitude. With the ratio of 
the coupling strength to the atomic frequency increas- 
ing, the times of entanglement revival decrease. There 
are entanglement revivals after a period of time of disen- 



tanglement, which are different from that of the case of 
(7$. 

(2) For ujQ — 27 and A = 2O7, the evolution dynamics 
of concurrence Cp is almost the same as that in Fig.8. 
Unlike the two cases above, the evolution behavior of 
concurrence Cq, becomes symmetric because the strong 
coupling of atom with non-Markovian reservoir and the 
effect of counter-rotating wave interaction. 

The above results reveal that, in principal, the decoher- 
ence of C$ is symmetrical with /3^ because of the sym- 
metry of initial state |$), while that is unsymmetrical 
with /3^ because the initial state |^) is unsymmetrical 
with /3^. However, the strong coupling and the counter- 
rotating wave interaction could make the decoherence of 
Ciji become symmetrical with /3^. With the enhancing 
of coupling strength, the memory effect of the counter- 
rotating wave terms becomes dominant because the the 
average value of correlation function a2 increases. The 
results are quite different from that of the case for qubits 



in a vacuum cavity [Iz 



IV. CONCLUSION 

The reduced non-perturbative quantum master equa- 
tion of atom in initially squeezed vacuum cavity has been 
derived and its exact solution is obtained. The decoher- 
ence behaviors of two qubits with squeeze parameters and 
other parameters have been discussed. 

The results show that the decoherence behavior of 
two qubits in a squeezed reservoir is dependent on the 
squeeze parameter, the ratio of the coupling strength 
to the atomic transition frequency and the ratio of the 
width of reservoir spectral density to the atomic transi- 
tion frequency. First, in principle, the disentanglement 
time decreases with the increasing of squeeze parame- 
ter r, due to the increasing of average photon number 



of every mode in squeezed vacuum cavity. Second, there 
are entanglement revivals after complete disentanglement 
for the case of even parity initial Bell state, while there 
are entanglement decrease and revival before complete 
disentanglement for the case of odd parity initial Bell 
state. Third, with the enhancing of coupling strength, 
the times of entanglement revival decrease due to the 
effect of counter-rotating wave terms. Fourth, there is 
entanglement revival due to the memory effect of the 
non-Mar kovian squeezed environment. 
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